Charged impurity scattering limited low temperature resistivity of low density silicon 

inversion layers 
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We calculate within the Boltzmann equation approach the charged impurity scattering limited low 
temperature electronic resistivity of low density n-type inversion layers in Si MOSFET structures. 
We find a rather sharp quantum to classical crossover in the transport behavior in the — 5K 
temperature range, with the low density, low temperature mobility showing a strikingly strong 
non-monotonic temperature dependence, which may qualitatively explain the recently observed 
anomalously strong temperature dependent resistivity in low-density, high-mobility MOSFETs. 
PACS Number : 73.40.-c, 73.40.Qv, 73.50.Bk, 71.30.+h 
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Several recent publications on low temperature resis- 
tivity measurements in various low density two 
dimensional (2D) systems report the observation of an 
anomalously strong temperature dependence as a func- 
tion of carrier density, which has been interpreted as 
evidence for a zero temperature two dimensional metal- 
insulator transition (2D M-I-T), which is considered to 
be forbidden in two dimensions (at least for a non- 
interacting 2D system) by the one parameter scaling the- 
ory of localization ||. A number of theoretical papers 
have appeared in the literature providing many pos- 
sible resolutions of this seemingly unanticipated (but ap- 
parently ubiquitous) phenomenon. In this Letter we pro- 
pose a possible theoretical explanation for (at least a part 
of ) the observed phenomena. Our explanation is quanti- 
tative, microscopic, and physically motivated. Although 
our theory is quite general and generic (and thus appli- 
cable to all the systems [Q-Q exhibiting the so-called 2D 
M-I-T), we specifically consider here the electron inver- 
sion layer in Si MOSFETs, which is both the original 
system in which the 2D M-I-T was first reported |jj and 
is also the most exhaustively experimentally studied 
system in this context. It is important to emphasize that, 
in contrast to much [[| of the existing theoretical work 
on the subject, our theory does not address the existence 
(or not) of a zero temperature 2D M-I-T, but specifically 
addresses the issue of quantitatively understanding the 
strikingly unusual finite temperature experimental results 
on the effective "metallic" side of the transition. 

We first summarize the key experimental features of 
the 2D M-I-T phenomenon (focusing on Si MOSFETs), 
emphasizing the specific aspects addressed in our theory. 
Experimentally one finds a "critical density" (n c ) sepa- 
rating an effective "metallic" behavior (for density n s > 
n c ) from an effective "insulating" behavior (n s < n c ). 
We concentrate entirely on the effective "metallic" behav- 
ior in this Letter since a 2D metal is "unusual" according 
to the conventional theory || and a 2D insulator is not. 
The experimental insulating behavior (for n s < n c ) is 



quite conventional for a strongly localized semiconductor 
and can be understood using standard transport models 
The effective "metallic" behavior is characterized 
by a strong drop in the temperature dependent resistiv- 
ity, p(T), at low temperatures (O.IK < T < 1 — 3K) 
and at low densities (n s > n c ). This novel and dramati- 
cally strong temperature dependence of p(T), where p(T) 
may drop by a factor of 2 — 10 at low electron densities 
as temperature decreases from 2K to 100 mK, is one of 
the most significant experimental observations we qual- 
itatively explain in this Letter. In addition the experi- 
mental resistivity, p(T, n s ), as a function of temperature 
and electron density shows an approximate "scaling" be- 
havior p(T,n s ) ~ p(T/Tq) with To = To(n s ) indicating 
consistency with quantum criticality. Our theoretical re- 
sults show the same "scaling" behavior with our calcu- 
lated Tq(ji s ) having very similar density dependence as 
the experimental observation. There are interesting as- 
pects of the magnetic field and the electric field depen- 
dence of the observed resistivity, which we do not address 
here, concentrating entirely on the behavior of p(T, n s ) 
in the n s > n c "metallic" regime. It is this "anomalous 
metallic" behavior (in the sense of a very strong metal- 
lic temperature dependence of the resistivity in a narrow 
density range above n c ) which has created the recent in- 
terest in the 2D M-I-T phenomena since in general, the 
temperature dependent resistivity of a metal should sat- 
urate as it enters the low temperature Bloch-Griineisen 
regime without manifesting any strong temperature de- 
pendence. 

Our theory, which provides good qualitative agree- 
ment with the existing experimental data on the metallic 
(n s > n c ) side of the transition, is based on two essential 
assumptions: (1) transport is dominated by charged im- 
purity scattering centers (with a density of Ni per unit 
area) which are randomly distributed at the interface; (2) 
the M-I-T at n s = n c is characterized by a "freeze-out" 
of free carriers due to impurity binding — the free car- 
rier density responsible for "metallic" transport is thus 
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(n s — n c ) for n s > n c , and on the insulating side, n s < n c , 
the free carrier density (at T = 0) is by definition zero. 
Some justifications for these assumptions have been pro- 
vided in ref. H although our current model transcends 
the specific scenario envisioned in ref. j7| and is more 
general. In contrast to ref. Q, we do not specify any 
particular mechanism for the carrier freeze-out and ac- 
cept it as an experimental fact. We note that we could 
extend our model and go beyond the above two assump- 
tions, for example, by making the effective free carrier 
density n = (n s — n c )9(n s — n c ) + n a (T), where n a (T) is 
a thermally activated contribution to the carrier density 
(this relaxes the second assumption), and/or by introduc- 
ing additional scattering mechanisms such as the short- 
range surface roughness scattering (this relaxes the first 
assumption). These extensions beyond our two essential 
approximations will undoubtedly produce better quan- 
titative agreement between our theory and experiment 
(at the price of having unknown adjustable parameters). 
We, however, refrain from such a generalized theory, be- 
cause we believe that the minimal theory, constrained by 
our two stringent assumptions and thus allowing for only 
one unknown parameter (the charged impurity density 
Ni) which sets the overall scale of resistivity in the sys- 
tem, already catches much of the essential physics in the 
problem. 

We use the finite temperature Drude-Boltzmann the- 
ory to calculate the ohmic resistivity of the inversion layer 
electrons taking only into account long range Coulom- 
bic scattering by the random static charged impurity 
centers with the electron-impurity Coulomb interaction 
being screened by the 2D electron gas in the random 
phase approximation (RPA). The resistivity is given by 
p = er -1 , where the conductivity a = ne 2 < t > /m with 
m as the carrier effective mass, and < r > is the energy 
averaged finite temperature scattering time: 



< T >- 



JdEEr(E) (-%) 
JdEE(-§L 



(1) 



where f(E) is the Fermi distribution function, f(E) = 
{1 + exp[(£' — /i)]/fcsT} _1 with p(T, n) as the finite tem- 
perature chemical potential of the free carrier system de- 
termined self-consistently. The energy dependent scat- 
tering time t(E) for our model of randomly distributed 
interfacial impurity charge centers is given by 
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with q = |k — k'|, 9 = 6>kk' is the scattering angle between 
k and k', E = e k = h 2 k 2 /2m, e k ' = h 2 k /2 /2m, v{q) is 
the 2D Coulomb interaction between an electron and an 
impurity, and e(q) = e(q; p,T) is the 2D finite tempera- 
ture static RPA dielectric (screening) function In 



calculating the Coulomb interaction and the RPA dielec- 
tric function in Eq. (2) we take into account subband 
quantization effects in the inversion layer through the 
lowest subband variational wavefunction ||. We note 
that there are two independent sources of temperature 
dependence in our calculated resistivity — one source is 
the energy averaging defined in Eq. (1) and the other is 
the explicit temperature dependence of the finite tem- 
perature dielectric function e(q) which produces a di- 
rect temperature dependence through screening in Eq. 
(2). At very high temperatures, when T 3> Tp with 
Tp = fi(T = 0)/feg as the free carrier Fermi temperature, 
the system is classical and it is easy to show that Eq. (1) 
leads to a decreasing resistivity with increasing temper- 
ature: p(T) ~ AiT/Tp)- 1 for T > Tp. In the quantum 
regime, however, energy averaging by itself produces a 
weak quadratic (negative) temperature dependence ac- 
cording to Eq. (1): p ~ p(T = 0) - B(T/T F ) 2 , for T < 
Tp. For Si inversion layer, however, this low temperature 
negative temperature dependence is overwhelmed |9|-[Tl|l 
by the temperature dependence of the screening function 
in Eq. (2) which dominates 2kp scattering — this phe- 
nomenon arises from the specific form of the 2D screening 
function which is a constant upto q = 2kp, and has a cusp 
at 2k f at T = 0. This strong temperature dependence 
arising from the low temperature screening function pro- 
duces a linear rise in the low temperature (T <C Tp) resis- 
tivity with increasing temperature according to Eq. (2): 
p(T) - p(T = 0) + C(T/T F ), for T -C Tp. This linear 
temperature dependence is, however, cut off at very low 
temperatures due to the rounding of the sharp corner in 
the 2D screening function by impurity scattering effects 
— at very low temperature T < T D where T D 
(= T/irkB with T as the collisional broadening) is the 
collisional broadening induced Dingle temperature, the 
explicit temperature dependence of e(q, T) is suppressed. 
At the densities and temperatures of interest in the 2D 
M-I-T phenomena all of these distinct physical effects are 
operational, and the actual behavior of p(T, n) could be 
quite complicated because the four different asymptotic 
mechanisms discussed above compete with each other as 
the system crosses over from a non-degenerate classical 
(T > Tp) to a strongly screened degenerate quantum 
(T << Tp) regime. We note that in general the tem- 
perature dependence is non-monotonic (particularly at 
lower densities where the energy averaging effects are 
significant), as has been experimentally observed 
because the temperature dependence of Eq. (1) by it- 
self produces a negative temperature coefficient whereas 
screening through Eq. (2) produces a positive tempera- 
ture coefficient. 

In Fig. 1 we show our numerically calculated resistiv- 
ity for the Si- 15 sample of ref. [JO. We show calculated 
p(T) as a function of T in Fig. 1 for several values of 
n s > n c . We use several different Dingle temperatures to 
incorporate ]10|-p2[ the impurity scattering induced 
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FIG. 1. The calculated resistivity p(T) for the Si- 15 sam- 
ple of ref. [1]. We use (a) the fixed Dingle temperature, 
Td = OK (dashed lines) and Tb = 0.5^ (solid lines), for var- 
ious electron densities, n s = 8.9, 9.2, 10.2, 12. 4x 10 10 cm~ 2 
(from top to bottom), and the experimental data points are 
taken from ref. [1]; (b) the density dependent Dingle tem- 
perature for various densities, n s = 8.9, 9.2, 10.2, 11.3, 11.8, 
12.4, 13.5, 14.5, 15.5, 16.5, 20.0, 24.2xl0 10 cm _2 (from top to 
bottom). 

collisional broadening corrections in the screening func- 
tion, including the pure RPA (Tb = 0) case which com- 
pletely neglects collisional broadening effects on screen- 
ing. In Fig. 1(b) we show the calculated Si- 15 results 
where the Dingle temperature varies as a function of elec- 
tron density. For each density the appropriate Tb (going 
into the screening calculation) is determined from the 
resistivity for that particular density. In Fig. 1(b) the 
temperature dependence of p(T) at low temperatures is 
strongest at intermediate densities somewhat away from 
n c whereas in Fig. 1(a) the temperature dependence of 
p(T) becomes stronger as one approaches n c , and is the 
strongest at the lowest density. This arises from the com- 
petition in screening among T, Tp, and Tb — at the 
lowest densities the temperature dependence is moder- 
ated by having relatively high values of Tb whereas at 
high densities the temperature dependence is suppressed 
by the large value of Tp, implying that the strongest 
temperature dependence of p(T) occurs at intermediate 
densities where neither Tb nor Tp is too high. Putting 
Tb = leads to stronger temperature dependence be- 
cause the temperature dependence of screening is not cut 
off at low temperatures as it is in the T < Tb < Tp 
regime for the Td ^ results. 

The impurity density has been fixed by demanding 



agreement between experiment and theory at high tem- 
peratures (T = 5K) and the highest densities. The im- 
purity density Ni thus sets the scale of the overall resis- 
tivity (p oc Ni), and does not affect the calculated T and 
n s dependence of p(T,n s ). It is important to emphasize 
that Ni values needed in our calculation to obtain quan- 
titative agreement with experiment are in the reasonable 
range of 2Vj ~ 10 10 cto~ 2 , which is known |?],[) 11 to be 
the typical effective random charged impurity scattering 
center density in high mobility Si MOSFETs. Since N 
is the only "free" parameter of our theory, it is signifi- 
cant that we obtain a reasonable value for Ni in order to 
achieve agreement between theory and experiment. We 
emphasize that our theory is valid even if the metallic 
(n s > n c ) system is actually weakly localized as long as 
the effective localization length is larger than the sys- 
tem size or the phase coherence length. In general, our 
calculated resistivity is higher (by 25 ~ 40%) than the 
experimental values at low densities (n s > n c ), i.e., our 
theory predicts a somewhat stronger n s dependence of 
p than that observed experimentally. This discrepancy 
can be corrected by adding an activated carrier density 
n a (T) to our effective carrier density n — (n s — n c ) +n a , 
which produces the strongest effect at the lowest densi- 
ties (and essentially no effect at higher densities), and 
would reduce p(T) at lower densities. One can also use 
a variable impurity density Ni(n s ) which varies with the 
gate voltage (following the spirit of ref. ||), and is lower 
at lower values of n s , again producing quantitative agree- 
ment between theory and experiment. Given the overall 
excellent qualitative agreement between our results and 
the experimental data of ref. , we think that these re- 
finements of our model are not particularly essential or 
meaningful. 

We show the experimental data points for Si- 15 taken 
from ref. [Q in Fig. 1(a) to give an idea about the level of 
agreement between our calculation and the experimental 
results. We do not attach particularly great significance 
to the quantitative agreement achieved in Fig. 1(a) be- 
cause of the various approximations in our theory. We 
do emphasize, however, that our calculations catch all 
the essential qualitative features of the low temperature 
experimental data We obtain the observed non- 

monotonicity in p(T) at low densities and also the strong 
drop in p(T) at low densities in the 0.1 ~ 2K temper- 
ature range. Consistent with the experimental observa- 
tions our calculated low density p(T) could drop by an 
order of magnitude for 1 — 2K change in the temperature. 
Our high density results show weak monotonic increasing 
p(T) with increasing T similar to experimental observa- 
tions We have carried out calculations for all the 
reported Si samples (as well as GaAs samples) in the lit- 
erature [Q-Q , and our level of qualitative agreement with 
experiment in uniformly good (typically as good as it is 
in Fig. 1) for all the existing experiments. For lower mo- 
bility samples, our calculated temperature dependence is 
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FIG. 2. The calculated scaling behavior of the resistivity 
from Fig. 1(b). Inset shows density dependence of the scaling 
parameter To. 

weaker (consistent with experimental findings) because 
screening is suppressed by stronger impurity scattering 
effects (through a higher value of level broadening or Din- 
gle temperature). 

In Fig. 2 we show our calculated "scaling" properties 
of p(T, n s ) ~ p(T/T ) with T = T (n s ) for the Si- 15 re- 
sults shown in Fig. 1(b). The scaling is obtained entirely 
numerically by obtaining the To which gives the best scal- 
ing fit to the calculated p(T,n s ). Comparing Fig. 2 with 
the corresponding experimental scaling plots Jl]] of resis- 
tivity we conclude that our theoretical scaling behavior 
of p(T, n s ) is about as good as the corresponding exper- 
imental scaling. In particular, our To(n s ), shown as an 
inset in Fig. 2, agrees reasonably well with the experi- 
mental results [Q. We obtain very similar "scaling" re- 
sults for the other Si samples in refs. The "scaling" 
we obtain in Fig. 2 underscores the important point that 
the experimentally observed scaling behavior in a narrow 
(T,n s ) range does not necessarily imply quantum criti- 
cality. 

Before concluding we point out the approximations 
made in our calculations. We have assumed uncritically 
that the Drude-Boltzmann transport theory, which is ex- 
tensively and successfully employed || in the device sim- 
ulation of Si MOSFETs, applies to the problem being 
studied. Our main justification for applying the stan- 
dard transport theory to the current problem is our belief 
that such a "zeroth order" , "one-parameter" (iVj being 
the only parameter in our model) theory must be applied 
to the problem and compared with the experimental data 
before one can discuss more speculative (and calculation- 
ally difficult) approaches |J . The fact that such a zeroth 
order theory already obtains good qualitative agreement 
with the experimental results indicates that charged im- 
purity scattering, carrier binding and freeze-out, temper- 
ature and density dependence of 2D screening, and clas- 
sical to quantum crossover (in the T = — 5K range) are 
playing significant roles in the experiments and cannot 



be neglected in any theoretical analysis of the "2D M-I- 
T" phenomenon. Our other approximations of using the 
RPA screening (we actually incorporate a 2D Hubbard 
local field correction ]T^ ] in our screening, which has no 
qualitative effect on our results) and the Dingle tempera- 
ture approximation to incorporate collisional broadening 
effects on screening are quite reasonable (at least qual- 
itatively) within our model and approximation scheme, 
and may be systematically improved (with a great deal of 
work) if future experiments warrant such a quantitative 
improvement of the theory. It is important to empha- 
size that quantum corrections, including localization ef- 
fects, are left out of our semi-classical Drude-Boltzmann 
theory. In providing some justification we mention that 
the effective dimensionless parameter kpl (where I is the 
transport mean free path) is typically 2 or larger in our 
results (at T — 0, which we believe to be the appropriate 
limit to consider), and therefore a Boltzmann theory may 
have reasonable validity. We estimate weak localization 
effects to be substantially weaker than the temperature 
dependence shown in Fig. 1 in the experimental temper- 
ature range (T > 50mK) in refs. jjjJI]. The fact that the 
observed temperature dependence, particularly at lower 
temperatures, is somewhat stronger than our calculated 
results may very well be the manifestation of quantum 
fluctuation or interaction effects neglected in our theory. 
An important approximation of our theory (consistent 
with the Drude-Boltzmann approach) is the neglect of 
inelastic electron-electron interaction, which may well be 
significant in the low density 2D systems of experimental 
relevance. For example, it is possible that the insulat- 
ing system (n s < n c ) is an electron glass (arising from 
the competition/frustration between interaction and dis- 
order). While a quantitative theory including disorder 
and interaction effects is extremely difficult, we speculate 
that our Boltzmann theory (in particular the quantum- 
classical crossover which leads to the strong temperature 
dependence) is sufficiently robust so that our qualitative 
conclusions will remain unaffected. 

We conclude by emphasizing the specific salient fea- 
tures of the qualitative agreement between our theory 
and experimental data on the metallic side: (1) strong 
temperature dependence at low and intermediate den- 
sities {n s > n c ); (2) non-monotonicity in p(T), aris- 
ing from quantum-classical crossover, at low values of 
n s > n c where p(T) increases weakly with decreasing T 
at higher temperatures and decreases strongly with T at 
lower temperatures; (3) scaling of p(T,n s ) ~ p(T/Tb) 
with the theoretical To(n s ) agreeing with the experi- 
mental results; (4) our calculated zero temperature con- 
ductivity, a(T — 0, n s ) = l/p(T — > 0, n s ), shows an 
approximately (within 25%) linear density dependence, 
cr(T = 0) oc n = (n s — n c ), which is consistent with the 
T — > extrapolation of the experimental M resistiv- 
ity and also with several other experimental [ 14j findings 
[this dependence, a{T = 0) oc (n s — n c ), also supports 
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our basic freeze-out or binding model]; (5) for increas- 
ing disorder (i.e., for lower mobility samples) we predict 
an increasing n c and weaker temperature dependence, as 
observed experimentally. 

This work is supported by the U.S.-ARO and the U.S.- 
ONR. 
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